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Endohedrally doped Si2o fullerenes appear as appealing building blocks for nanoscale materials. 
We investigate their structural stability with an unbiased and systematic global geometry opti- 
mization method within density-functional theory. For a wide range of metal doping atoms, it was 
sufficient to explore the Born Oppenheimer surface for only a moderate number of local minima 
to find structures that clearly differ from the initial endohedral cages, but are considerably more 
favorable in terms of energy. Previously proposed structures are thus all met ast able. 
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As miniaturization techniques are reaching their ulti- 
mate limits, the interest in novel silicon based nanoscale 
devices increases. Notwithstanding, most common ma- 
terials in nano sciences to date are carbon based such as 
fullerenes and nanotubes. Si based clusters and nanopar- 
ticles have also been studied extensively and it was shown 
that their electronic structure as well as their mechani- 
cal, optical and magnetic properties can be manipulated 
by changing their shape, size and composition. There is 
widespread hope that such Si nanomaterials may be basic 
building blocks for more complicated structures, such as 
wires and layers [l|-[5| . Unfortunately, up to now no stable 
Si building blocks have been found that are as chemically 
unreactive and symmetric, and therewith attractive for 
cluster assembled materials, as the carbon fullerenes. 

Endohedral doping with metal atoms is a primary av- 
enue believed to stabilize cage-like Si geometries. In 
fact, clathrates are composed of corresponding polyhe- 
dral building blocks [6]. The exceptional elastic, ther- 
moelectric, optoelectronic and super-conducting proper- 
ties of these porous crystals [7-10] already illustrate the 
unique potential offered if novel materials could be tai- 
lored out of such Si-based subunits. Considering that 
C20 forms the smallest known fullerene, Si2o clusters 
represent a particularly interesting size in this context 
that should in principle be large enough to encapsulate 
a metal atom [11]. In contrast to the intrinsically unsta- 
ble hollow Si2o fullerene [12], endohedral doping with a 
range of metal atoms was indeed theoretically predicted 
to stabilize the cage structure [l3|, [l^ • 

In these, as well as in numerous equivalent theoretical 
studies on other cluster sizes, the stability was inferred 
from computed embedding and binding energies of re- 
laxed structures. By construction, corresponding geom- 
etry optimizations lead, however, only to the next lo- 
cal minimum on the Born- Oppenheimer potential energy 
surface (PES). While a harmonic frequency analysis may 



ensure that this local minimum has indeed been reached, 
this still does not tell anything about the global PES fea- 
tures. In particular such an approach does not tell us 
if there are other energetically even more favorable min- 
ima, or if the present structure indeed corresponds to 
the global minimum. Starting the geometry optimization 
from several initial configurations [15, 16] or using several 
stages of symmetry constraints 0,ll8| may provide some 
information in this direction. Still, the corresponding ex- 
ploration of the PES is by no means systematic, and the 
rehability of the deduced structural stability uncertain. 

In this work, we therefore reexamine the structure of 
metal-doped Si2o clusters using a global and unbiased ge- 
ometry optimization technique within density-functional 
theory (DFT). For essentially the entire range of previ- 
ously proposed metal dopants this readily identifies sig- 
nificantly more stable structures that no longer corre- 
spond to endohedral fullerene cages. The latter config- 
uration thus only corresponds to a local PES minimum, 
and the partial information we obtain on the surround- 
ing barriers even suggests that this minimum is in most 
cases quite shallow. With a corresponding at best fee- 
ble metastability restricted to low temperatures, doped 
Si2o fullerenes are unlikely useful building blocks for fu- 
ture nanoscale materials - unless additional stabilization 
mechanisms are identified. 

The minima hopping method (MHM) [l9|] is designed 
to find the global minimum of complex polyatomic sys- 
tems in an efficient way. The general idea is to limit re- 
peated visits of the same local minima without penalizing 
crossings through important transition regions, such as 
hubs connecting superbasins of the potential landscape. 
The method is composed of an inner part, that attempts 
to escape from the current minimum by following short 
trajectories from molecular dynamics (MD), and an outer 
part that either accepts or rejects the new configuration 
by simple energy thresholding. A feedback mechanism 
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TABLE L Results of the MHM search for the lowest-energy 
structure of doped Si2o clusters, which either correspond to 
exohedral configurations with broken cages or smaller Sin 
cages with peripheral bud. Additionally shown is the en- 
ergy difference (in eV) between this most stable structure en- 
countered and the lowest-energy intact cage (see text). Sum- 
marized is the energetic data for a range of local, gradient- 
corrected and hybrid DFT xc functionals as obtained with 
FHLaims [28]. The listed energies are always the energies 
that correspond to the lowest energy spin state. 



on both parts allows to take advantage of the history of 
minima visited, as well as of the Bell-Evans-Polanyi prin- 
ciple, which correlates lower energy barriers with deeper 
basins (2q| . 

In order to reach predictive quality, the PES explored 
by the MHM must be computed from first-principles. 
Here we use DFT as implemented in the BigDFT pack- 
age [2H with valence-type pseudopotentials for this 
purpose. In the spirit of the dual MHM [23| two lev- 
els of accuracy are considered to reduce the computa- 
tional cost. During the MD escapes and in the initial 
stages of local geometry relaxations a coarser grid with 
smaller simulation boxes was chosen to define the em- 
ployed adaptive wavelet basis. For the final geometry 
optimization and the evaluation of the total energy of 
the relaxed structure highly accurate parameter sets were 
used. In these calculations we rely on the widespread 
local-density approximation (LDA) as an efficient 
general-purpose approach to treat electronic exchange 
and correlation (xc). In order to check the accuracy of 
the LDA xc functional we recomputed the energetic order 
of the identified minima with gradient-corrected (PBE 
^]) and hybrid (PBEO [26,], B3LYP \^]) functionals. 
The latter computations were done with the accurate 
all-electron full-potential code FHI-aims [28] using the 
"tier2" basis-set composed of atomic-centered numeric 
orbit als. For the LDA and PBE functionals contained in 
both codes the obtained energetic differences agreed to 
within 150 me V, thereby confirming the accuracy of the 
pseudopotentials employed in the initial BigDFT calcu- 
lations and the near-completeness of the basis set used 
in the FHI-aims calculations. 



A natural starting point for our investigation of the 
structural stability of metal-doped Si2o fullerenes are the 
perfect cages of Ih symmetry as illustrated in Fig. [TJ 
As impurity atoms we consider a wide range of sim- 
ple and transition metals (Ba, Ca, Cr, Cu, K, Na, Pb, 
Kb, Sr, Ti, V and Zr), which comprises in particular 
those elements contained in clathrate materials and those 
previously proposed to stabilize the Si2o cage structure 
|l3l,[l^. Surprisingly, only relatively short MHM searches 
over a few hundred PES minima were necessary for all 
of these dopants to reveal significantly more stable struc- 
tures that deviate in either of two ways qualitatively from 
an endohedral fullerene configuration: As summarized in 
Fig. [1] and Table HI for most impurity atoms exohedral 
structures were readily identified. In all other cases, the 
dopant was not expulsed, but encapsulated in a smaller 
cage with the remaining Si atoms forming an apical bud. 
For Ti, V, Cr and Cu dopants, this endohedral struc- 
ture of the lowest-energy configuration is fully consistent 
with the interpretation of Ar physisorption experiments 
ll|. The specific size of the identified smaller cages is 



furthermore in line with a preferred stability of the cor- 
responding CuSiio, CrSii5, and MSiie (M =Ti, V, Zr) 
clusters as deduced from their abundance in mass spec- 
tra or simple electron counting rules [ll|, [29|, [30|. TiSiie 
and TiSii6 cages have also already been identified as local 
minima in DFT calculations (3^ . 

The MHM searches were stopped as soon as configu- 
rations of significantly lower energy than the initial sym- 
metric Ih cage were identified. During the runs typically 
also a number of more favorable configurations were vis- 
ited, in which the cage was (partly heavily) distorted, 
but could still be considered intact, cf. Fig. [TJ For those 
dopant atoms ultimately leading to completely broken 
cages, a distinction of these structures from the lowest- 
energy exohedral ones is rather unambiguous. As shown 
in Table H] the corresponding energy gap to the most sta- 
ble of these identified intact cage structures is in all cases 
quite large. This holds for all of the employed xc func- 
tionals, even though quite some quantitative scatter can 
be discerned. With such a clear gap, it is unlikely that 
an intact cage structure would exist that is even lower 
in energy and that has been missed in the performed fi- 
nite MHM searches. Instead, we rather expect that in 
analogy to pure Si clusters there exists a multitude of 
further disordered exohedral configurations, which are all 
extremely close in energy to the here identified most fa- 
vorable structure [3l|. With exohedral cages being also 
favorable in terms of entropy, we therefore conclude that 
the hitherto proposed endohedral Si2o fullerenes for the 
corresponding dopant atoms are only met ast able. 

For those impurity atoms ultimately encapsulated by 
a smaller number of Si atoms, already a mere relaxation 
of the initial Ih cage resulted in rather heavy distortions 
as illustrated in Fig. [TJ Subsequently sampled configu- 
rations exhibited more and more pronounced distortions. 
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FIG. 1: The energetically lowest-lying configurations found for all M@Si2o systems as shown in the second and fourth columns 
have either lost their cage structure and dopant (bright) or undergone a transition to a smaller distorted cage with some 
peripheral silicon atoms. The corresponding most stable (distorted) cage-like structures are given in columns one and three. 
From top to bottom: M — Ba, Ca, Cr, Cu, K, and Na on the left, and M — Pb, Rb, Sr, Ti, V and Zr on the right half. 
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spanning a rather continuous range up to the smaher cage 
lowest-energy structure. In this situation the specifica- 
tion of an energy gap to the lowest-energy intact Si2o cage 
is not well defined and we therefore quote in Table |I] the 
energy difference to the initial relaxed Ih cage. Again, 
this energy difference is sizable in all employed xc func- 
tionals. From this and the observed range of increasingly 
distorted cages we would therefore also rule out for these 
dopants that more stable fullerene configurations exist 
that were not identified in the present MHM searches. 

The thus disclosed metastability of the endohedral Si2o 
cages for a wide range of dopant atoms is in distinct con- 
trast to the high-symmetry carbon fullerenes [32j. The 
latter correspond to the global PES minimum with a 
large energy gap to the next lowest-energy structures 
formed by point defects. Still, if the barriers surround- 
ing the local PES minima corresponding to the endohe- 
dral cages are sufficiently high, kinetic trapping might 
be sufficiently long. However, it is well known that clus- 
ters of low symmetry have a broad distribution of barrier 
heights [33I and that dynamic processes in such systems 
involve in general the crossing of several barriers. It is 
therfore unlikely that low symmetry cage minima are sur- 
rounded only by high barriers. Exploring a dozen saddle 
points of the CaSi2o cage we found barriers that were in 
most cases higher than 1 eV. One barrier leading to an 
opening of the cage was even only 0.85 eV. According to 
kinetic rate theory such a barrier would correspond to a 
life time not longer than a few seconds at room tempera- 
ture. Of course, this showcase does not allow to exclude 
kinetic trapping for all studied dopants in general. It nei- 
ther provides a complete pathway from the cage to the 
exohedral structure. Even though more detailed studies 
of the dynamics of these clusters would be required to 
precisely predict their life time, our present results sug- 
gest that larger barriers than those identified here would 
be required to stabilize the metastable CaSi2o cage over 
time spans relevant for materials applications. 

In conclusion, we have used DFT based global geome- 
try optimization to reexamine the proposed stabilization 
of Si2o fullerenes through endohedral metal doping. For 
a wide range of simple and transition metal dopants this 
readily reveals that the desirable symmetric cage struc- 
ture is only metastable. Either exohedral compact config- 
urations or endohedral smaller cages with excess Si atoms 
forming an apical bud are instead significantly more sta- 
ble. Regardless of whether local, gradient-corrected or 
hybrid DFT xc functionals are employed, the resulting 
energy gap of these lowest-energy configurations to the 
metastable fullerene cage is in most cases in excess of 
leV. 

These findings put severe doubts on the dream of sili- 
con based fullerenes as the building blocks for nano sci- 
ences. Beyond the specific Si2o cages examined here, 
our study furthermore underscores the importance of a 
systematic exploration of the configurational space when 
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